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ABSTRACT 

Stiffened  plates  are  designed  for  minimum  volume  subject 
to  constraints  on:   Von  Mises  maximum  stresses,  nodal  dis- 
placements, height  to  thickness  ratios  of  frame  webs,  and 
width  to  thickness  ratios  of  frame  flanges.   Design  variables 
are  plate  thicknesses  and  stiff ener  dimensions. 

A  finite  element  analysis  program  is  developed  for  the 
design  of  stiffened  plates  using  numerical  optimization 
techniques.   The  program  may  be  used  as  a  stand  alone  analy- 
sis tool  or  may  be  coupled  to  an   optimizer  of  user's  choice. 
Rectangular  plate  elements  and  frame  elements  are  used  for 
the  idealization  of  stiffened  plates. 

Design  examples  are  presented  tC'  demonstrate  the  design 
method. 
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I.   INTRODUCTION 

The  development  of  high-:3peed  digital  computers  has  made 
possible  significant  changes  in  the  structure  design  process. 
One  of  these  is  the  availability  of  various  mathematical 
programming  methods  for  use  in  design  optimization.   The 
computer's  speed  allows  the  designer  to  now  consider  a  much 
wider  range  of  design  alternatives.   The  optimization  proce- 
dure provides  a  means  of  systematically  choosing  from  among 
these  alternatives  based  on  some  predetermined  rational 
criterion. 

Even  when  the  selected  numerical  method  is  able  to  arrive 
at  the  optimum  design^  tne  rssult  is  only  as  good  as  the 
design  model.   Here  there  is  an  even  greater  need  for  design 
experience  and  sound  engineering  judgment.   The  design  model 
must  be  carefully  developed  to  realistically  represent  the 
design  in  question. 

The  finite  element  model  for  stiffened  plate  is  shown  in 
Fig.  1.1,  where  the  plate  has  been  descritized  by  rectangular 
plate  elements  and  the  stiff ener  by  frame  type  elements.   The 
eccentricity  of  stiffener  is  transformed  to  the  linked  nodal 
point  by  applying  a  linear  equation  that  realtes  displacement 
degrees  of  freedom. 

The  purpose  of  this  thesis  is  to  develop  a  finite  element 
analysis  program  for  stiffened  plates,  and  to  design  the 
optimum  stiffened  plate  by  coupling  two  programs;  the 
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Figure  1.1    Stiffened  Plate  Model 
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analysis  program  andageneral  purpose  non-linear  optiimizer , 
COPES/CONMIN  [Refs.  1-2]. 

The  remaining  chapters  of  this  thesis  are  outlined  as 
follows : 

Chapter  II  presents  briefly  the  finite  element  method 
used  in  the  analysis  program  for  stiffened  plates. 

Chapter  III  presents  the  basic  concepts  of  the  cptin.i- 
zation  methods  used  in  the  COPES/CONMIN. 

Chapter  IV  presents  design  examples. 

Chapter  V  offers  conclusions  and  recommendations. 
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II.       FINITE    ELEJyiENT   METHOD 

A.       GENERAL 

The  finite  element  method  is  now  well  established  as  an 
engineering  tool  of  wide  application.   The  fundamental 
concept  of  the  finite  element  method  is  that  any  continuous 
quantity,  such  as  pressure  or  displacements,  can  be  approxi- 
mated by  a  discrete  model  composed  of  a  set  of  piecewise 
continuous  functions  defined  o-^rer   a  finite  number  of  sub- 
domains.   The  piecewise  continuous  functions  are  defined 
using  the  values  of  the  continuous  quantity  at  a  finite 
number  of  points  in  its  domain . 

The  formulation  of  the    fin.te  element  method  can  be  traced 
to  energy  procedures,  principally  the  minimum  potential  energy 
method  and  the  minimum  complimisntary  energy  method.   The 
minimum  potential  energy  method  is  associated  with  assumed 
displacement  parameters  as  unk]iowns  and  is  usually  termed 
the  "displacement"  or  "stiffness"  method.   On  the  other  hand, 
the  minimum  complimentary  energy  method  dealt  with  stress 
parameters  and  is  termed  the  "flexibility"  or  "force"  approach 
The  ease  with  which  a  continuous  displacement  pattern  can  be 
prescribed  (compared  to  the  alternative  approach  of  forming 
an  equilibrating  internal  force  field)  has  aided  the  wide- 
spread use  and  development  of  the  finite  element  displace- 
ment approach.   The  displacement  model  and  the  stiffness 
analysis  are  employed  in  the  analysis  program  developed  here. 
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This  chapter  will  briefly  present  some  of  the  general 
concepts  of  the  finite  element  method  used  in  the  analysis 
program. 

B.   FINITE  ELEMENT  DISPLACEMENT  APPROACH 

The  displacement  formulation  involves  derivation  of  the 
stiffness  matrix  of  each  individual  element.   The  stiffness 
matrix  of  the  entire  assembled  structure  is  then  obtained  by 
super-position.   This  matrix,  along  with  the  prescribed 
displacement  boundary  conditions  and  loads,  is  used  for  the 
solution  of  displacements  and  stresses. 

1.   Element  Analysis 

For  the  structural  applications  at  least,  the 
governing  equilibrium  equations  can  be  obtained  by  mini- 
mizing the  total  potential  energy  of  the  system.   The  total 
potential  energy,  it  can  be  expressed  as 

TT   =   y  /   a"^  £  dV  -  /   u*^  f  dV  -  /   u*^  q  dS      (2.1) 

where  a   and  e  are  the  stress  and  strain  vectors  respectively, 
u  the  displacements  at  any  point,  f  the  body  forces  per  unit 
volume  and  q  the  applied  surface  tractions.   Integrations  are 
taken  over  the  volume  V  of  the  structure  and  loaded  surface 
area  S. 

The  first  term  on  the  right  hand  side  of  Eq.  (2.1) 
represents  the  internal  strain  energy  and  the  second  and 
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third  terms  are  respectively  the  work  contributions  of  the 
body  forces  and  distributed  surface  loads. 

In  the  finite  element  displacement  method,  the  basi« 
steps  for  derivation  of  the  element  stiffness  matrix  are: 
a.   Express  the  displacements  to  have  unknown  values  only 
at  the  nodal  points,  so  that  the  variation  within  any 
element  is  described  in  terms  of  the  nodal  values  by 
means  of  interpolation  functions.   Thus 


u   =   N  u®  (2.2) 


where  N  is  the  set  of  interpolation  functions  termed 
the  shape  functions  and  u   is  the  vector  of  nodal 
displacements  of  the  element, 
b.   Express  the  strains  within  the  element  from  the 
element  nodal  displacements  as 


£   =   B  u®  (2.3) 


where  B  is  the  strain-displacement  matrix  generally 
composed  of  derivatives  of  the  shape  functions. 

c.  Express  the  stresses  relating  to  the  strains  by  use  of 
an  elasticity  matrix  D  containing  the  appropriate 
material  properties,  as  follows: 

a   =   D  £  (2.4) 

d.  Establish  the  equilibrium  equation  of  element. 
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Provided  that  the  element  shape  functions  have  been  chosen 
so  that  no  singularities  exist  in  the  integrands  of  the 
functional,  the  total  potential  energy  of  the  continuum  will 
be  the  sum  of  the  energy  contributions  of  the  individual 
elements.   Thus 


TT    =    y  TT  (2  .  d) 

^   e 
e 


where  71   represents  the  total  potential  energy  of  element  e 
which,  on  use  of  Eq.  (2.1),  can  be  written 


TT    =   t  /    u^*^  B*^  D  B  u^  dV  -  /    u^*^  N*^  f  dV 
e  e 

-  /    u^*^'  N^  q  dS  (2.6) 

S   ~    z      -         ^ 
e 

where  V   is  the  element  volume  and  S   the  loaded  element 
e  e 

surface  area.   Performance  of  the  minimization  for  the 
element  e  with  respect  to  the  noc.al  displacements  u^  f  or  the 
element  results  in 


=   j     (B  DB)u  dV   -  /    N  f  dV  -  /    N  q  dS 

3u        v^~~       ^     V=S~ 
e  e  e 


k^  u^  -  p^  (2.7) 


where 
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p®   =   /    n'^  f  dV   +  /    N*^  q  dS  (2.8) 


V   -   ~    e     S   -   -    ^ 
e  e 


are  the  equivalent  nodal  forces  for  the  element,  and 


k    =   /    B   D  B  dV  (2.9) 


is  termed  the  element  stiffness  matrix. 
2 .   Direct  Stiffness  Method 

The  real  elastic  structure  is  now  represented  by  a 
finite  number  of  small,  discrete  elements.   Once  their 
approximate  behaviors,  identified  by  their  individual  stiff- 
ness matrices  k^of  Eq.  (2.9),  have  been  established,  the 
sticfness  matrix  K  for  the  complete  structure  is  obtained  by 
the  proper  summation  of  each  element  stiffness  matrix  in  the 
structure.   The  summation  of  the  terms  in  Eq.  (2.7)  over  all 
the  elements,  when  equated  to  zero,  results  in  a  system  of 
equilibrium  equations  for  the  complete  continuum.   This 
assembly  process  is  known  as  the  Direct  Stiffness  Method. 
These  equations  are  then  solved  by  any  standard  technique  to 
yield  the  nodal  displacements.   Note  that  K  is  symmetric  and 
positive-definite. 

C.   FINITE  ELEMENTS  USED 

The  finite  elements  used  in  the  analysis  program  will  be 
described  briefly  in  this  section. 
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1.  Rectangular  Plate  Element 

The  rectangular  plate  element  used  here  is  illus- 
trated in  Fig.  2.1,  where  each  nodal  point  has  6  degrees  of 
freedom,  3  transitional  displacements  and  3  rotational 
displacements.   This  element  with  4  corner  nodal  points  has 
the  element  stiffness  matrix  of  order  24  x  24.   The  corres- 
ponding displacements  of  each  node  u,  , u^ , . . . , u^ .  will  be 
taken  to  be  positive  in  positive  directions  of  the  xyz- 
coordinates . 

The  main  assumptions  in  the  method  are  that  displace- 
ments are  small  compared  with  plate  thickness,  the  stress 
normal  to  the  midsurface  of  the  plate  is  negligible,  and 
normals  to  the  midsurface  before  deformation  remain  straight 
but  not  necessarily  normal  to  the  midsurface  after  deformation 

The  assumed  displacement  functions  will  be  taken  to 
be  linearly  varying  in  the  plane  of  the  element.   These 
displacement  functions  will  ensure  both  deflection  and  slope 
compatibility  of  the  adjacent  elements.   The  stiffness  matrix 

of  the  rectangular  plate  element,  which  is  matrix  product 

T 
B  DB  of  Eq.  (2.9)  integrated  over  its  volume,  is  summarizea 

by  Przemieniecki  [Ref .  3]  . 

Note  that  element  stiffness  matrices  are  formed 
directly  in  the  global  coordinate  system  so  that  no  trans- 
formations from  local  to  global  coordinates  are  required. 

2 .  Frame  Element 

The  frame  element  as  a  stiffener  has  6  degrees  of 
freedom  for  each  node,  such  as  those  of  plate  elements.   The 
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Figure  2.1    Rectangular  Plate  Element 
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basic  assumptions  are  that  the  stiffener  element  is  straight 
with  uniform  cross  section  capable  of  resisting  axial  forces, 
bending  moments  about  the  two  principal  axes  in  the  plane  of 
its  cross-section,  and  twisting  moments  about  centroidal  axis, 
and  that  its  deflection  due  to  shearing  strains  are  neglected. 
The  stiffener  element  has  2  kinds  of  options  which  are  illus- 
trated in  Fig.  2.3.   The  width  and  thickness  of  flange  of 
the  rectangular  stiffener  element  may  be  referred  as  zeros. 

In  order  to  determine  the  stiffness  property  of  a 
complete  structure,  a  global  coordinate  must  be  established 
for  all  unassembled  structural  elements  so  that  all  the 
displacements  and  their  corresponding  forces  will  be  referred 
to  this  system.   Since  the  element  stiffness  matrices  k  are 
initially  calculated  in  local  coordinates,  suitably  oriented 
to  minimize  the  computing  effort,  it  is  necessary  to  intro- 
duce transformation  matrices  changing  the  frame  of  reference 
from  a  local  to  a  global  coordinate  system.   The  first  step 
in  deriving  such  a  transformation  is  to  obtain  a  matrix  rela- 
tionship between  the  element  displacements  u  in  the  local 
system  and  the  element  displacements  u'  in  the  global  system. 
This  relationship  is  expressed  by  the  matrix  equation 

u  =   T  u'  (2.10) 

where  T  is  a  matrix  of  coefficients  obtained  by  resolving 
global  displacements  in  the  directions  of  local  coordinates. 
The  transformation  matrix  T  is  given  by 
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Figure  2.2    Frame  Element 
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a.  Rectangular  type  fr?jne  stiffener 


b.  T  type  frame  stiffener 


Figure  2.3    Stiffener  Options 
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repres.ents  matrices  of  direction  cosines  for  local  x,  y  and 
z  directions,  respectively,  measured  in  global  system  x,  y 
and  z.   Similarly  the  relationship  for  forces  is  expressed 
by  the:  matrix  equation 


p   =   T  p' 


(2.11) 


where  p  is  the  force  vector  in  the  local  system  and  p'  is 

the  force  vector  in  the  global  system.   Matrix  A  is  orthogonal; 

-1     T 
that  is,  T   =  T  .   Therefore 


u 


=   T   u 


(2.12) 


and 


p'   =   T^p 


(2.13) 
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Let  u*  and  u  be  two  ways  to  describe  the  same  virtual  dis- 
placement.  Virtual  work  is 

5u''^p'   =   6u'^p   =   6u''^t'^p  (2.14) 


So 


6u''^(p'  -  t'^P)   =   0,     P'   =   t'^P  (2.15) 


To  transform  the  stiffness  matrix,  we  start  with  k'u'  =  p' 
and  substitute  from  the  preceding  equations. 

k'u'    =   p'    =   t'^p    =   T'^k  u   =   T'^k  Tu'  (2.16) 

Since  this  relation  is  presumed  valid  for  any  u',  we  conclude 
that  the  required  stiffness  transformation  is 


T*^  k  T  (2.17) 


D.   OFFSETTING  OF  RIGID  LINKS 

One  of  the  most  important  advantages  of  the  finite  ele- 
ment technique  is  that  an  assembly  of  different  structure 
elements  such  as  plates  and  frames  can  be  dealt  within  a 
single  coordinate  analysis.   Usually  the  neutral  siarfaces  of 
plate  and  stiff ener  are  not  coincident:   the  stiff ener  is  on 
one  side  of  the  plate.   A  standard  preliminary  treatment  is 
to  connect  adjacent  plate  and  stiffener  nodes  by  a  rigid 
link,  so  that  degrees  of  freedom  of  the  stiffener  are  replaced 
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by  degrees  of  freedom  of  the  plate.   The  usual  assembly  is 
then  possible.   The  necessary  transformation  is  now  described. 
The  stiff ener  element  in  Fig.  2.4  has  the  usual  12  degrees 
of  freedom--6  at  node  A  and  6  at  node  E. .   With  reference  to 
these  degrees  of  freedom,  element  load  and  stiffness  matrices 
are  p  and  k.   Similar  degrees  of  freedom  are  used  at  nodes  1 
and  2  of  plate  element  of  rigid  links  A-1  and  B-2.   The 
"master"  degrees  of  freedom  at  node  1  and  "slave"  degrees 
of  freedom  at  node  A  have  the  relation 


f   u 


Al 


u 


A2 


u 


A3 


u 


A4 


u 


A  5 


u 


A6  J 


=   A. 


u 
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u 


u 
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u 


15 


u 
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(2.18) 


where 
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plate  element 


frame  element 


Figure  2.4    Rigid  Links 
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A  similar  expression  is  written  for  link  B-2  by  replacing 
subscripts  A  and  1  by  B  and  2.   The  transformation  vectors 
p  and  matrix  k,  associated  with  degrees  of  freedom  at  nodes  1 
and  2 ,    are 


p   =   T  p 


(2.19) 


k   =   T*^  k  T 


(2.20) 


where 


:12xl2 


^1 


0  ^ 


^2. 


The  foregoing  transformation  makes  the  translational 
displacements  depend  on  the  rotational  displacements  and  so 
introduces  an  unwanted  quadratic  field  into  the  translational 
displacements.   It  is  found  that  a  frame-stiffened  plate  is 
overly  flexible  in  a  coarse  mesh.   Mesh  refinement  helps: 
error  is  reduced  by  a  factor  of  four  if  the  number  of  elements 
is  doubled.   The  error  can  be  eliminated  through  the  addition 
of  one  more  degree  of  freedom.   Further  details  of  the  error 
reduction  are  described  in  [Ref.  4]. 
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III.   OPTIMIZATION  TECHNIQUES 

A.  INTRODUCTION 

In  this  chapter  some  of  the  fundamental  ideas  and  formu- 
lation methods  of  mathematical  progranmiing  are  introduced 
to  understand  the  optimization  method  and  the  optimizer 
COPES/CONMIN.   Design  parameters  used  here  are  then  dis- 
cussed.  Fox  [Ref.  5]  and  Himmelblau  [Ref .  6]  provide  an 
extensive  discussion  of  numerical  optimization  techniques 
and  their  application  to  engineering  design. 

B.  DEFINITIONS 

In  discussing  the  optimization  methods,  the  following 
definitions  will  be  useful: 

Design  variables — the  design  variables  are  the  numerical 
parameters  for  which  values  are  to  be  chosen  in  producing  a 
design.   In  a  structural  problem,  they  might  be  plate  thick- 
ness, frame  dimensions,  etc. 

Objective  function — the  objective  function,  is  the  single 
valued  function  with  respect  to  which  the  design  is  optimized. 
In  a  structural  design  problem,  it  might  be  the  weight,  volume 
or  fabricated  cost  of  the  structure.   The  selection  of  an 
objective  function  can  be  one  of  the  most  important  decisions 
in  the  optimum  design  process. 

Constraints — practical  design  problems  are  usually  sub- 
ject to  a  series  of  constraints  which  must  be  satisfied  in 
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order  to  produce  an  acceptable  design.   These  constraints 
may  be  linear  or  non-linear.   If  a  parameter  is  beyond  the 
value  of  a  specified  value,  it  is  said  to  be  violated. 

Side  constraint — the  side  constraint  is  a  constraint 
which  restricts  the  specified  range  of  a  design  variable 
for  reasons  other  than  the  direct  consideration  of  performance 

Feasible  design — the  feasible  design  (or  acceptable 
design)  is  a  design  which  satisfies  the  specified  constraints. 

Infeasible  design — the  infeasible  design  (or  unaacep- 
table  design)  is  a  design  which  does  not  satisfy  the  speci- 
fied constraints. 

The  general  nonlinear  optimization  problem  is  expressed 
mathematically  as  follows: 

Minimize  F(x)  (3.1) 

Subject  to: 


g.(x)  <_     0    for   j  =  l,m  (3.2) 


X.   <   X.   <   X.    for   i  =  l,n  (3.3) 

1—1—1 


where  the  vector,  x,  is  the  vector  of  n  design  variables. 
The  objective  function,  F(x),  given  by  Eq.  (3.1),  as  well 
as  the  constraint  functions  given  by  Eq.  (3.2),  may  be  linear 
or  non-linear  functions  of  the  design  variables.   They  may  be 
explicit  or  implicit  functions  of  x,  but  must  have  continuous 
first  derivatives.   If  it  is  desired  to  maximize  F(x),  the 
minimization  of  the  objective  function  is  used  since  maximum 
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F(x)  can  be  treated  miniinuin  of  -F(x)  .   The  function  g.  (x) 
is  the  set  of  inequality  constraints  to  be  met.   The  m  in- 
equality constraints  must  be  satisfied  to  be  a  feasible 

I  u 

design.   Side  constraints,  x.  and  x.,  are  the  lower  and  upper 

limitations  placed  on  the  design  variables.   Side  constraints 

could  be  included  in  Eq.  (3.2),  but  are  treated  separately 

for  efficiency. 

The  n-dimensional  space  spanned  by  the  design  variables 
X.  is  referred  to  as  the  design  space.   As  stated  previously, 
any  design  which  satisfies  the  inequalities  of  Eq.  (3.2)  is 
referred  to  as  a  feasible  design.   If  the  design  violates  one 
more  of  the  inequalities,  it  is  said  to  be  infeasible.   The 
minimum  feasible  design  is  said  to  be  optimal. 

Most  nonlinear  optimization  programs  update  the  vector 
of  design  variables  by  the  iterative  relationship 

X^+l   =   X^  +  a*  S^  (3.4) 

where  q  is  the  iterative  number,  vector  S  is  the  direction 
of  search  in  the  design  space,  and  a*  is  a  scalar  which 
defines  the  distance  of  travel  in  the  direction  S  during  the 
qth  iteration.   An  initial  design  defined  by  X  must  be 
defined  and  may  be  a  feasible  or  infeasible  design. 

The  optimization  process  then  proceeds  in  two  steps.   The 
first  is  the  finding  of  S  which  will  improve  the  design 
without  violating  constraints  and  the  second  is  the  determina- 
tion of  a*  which  will  improve  the  design  as  much  as  possible 
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when  moving  in  this  direction.   The  process  is  then  repeated 
until  there  is  no  further  design  improvement,  indicating 
that  this  is  the  optimum  attainable  design. 

C.   COPES/CONMIN 

The  COPES/CONMIN  optimization  progrcim  is  a  general  pur- 
pose, non-linear  optimizer  capable  of  handling  large,  con- 
strained problems.   This  includes  the  conjugate  direction 
method  of  Fletcher  and  Reeves  [Ref.  7]  for  unconstrained 
function  minimization  and  a  modification  of  Zoutendijk's 
Method  [Ref.  8]  of  Feasible  Directions  for  constrained  func- 
tion minimization.   It  has  been  successfully  used  in  connec- 
tion with  structural  optimization  [Ref.  9],  airfoil  design 
[Ref.  10],  aircraft  synthesis  [Ref.  11],  and  numerous  other 
engineering  applications. 

It  was  necessary  to  develop  a  subroutine,  ANALIZ,  which 
for  a  given  design,  would  analyze  a  stiffened  plate,  and  which 
would  be  suitable  for  coupling  with  the  optimizer.   The 
common  block  GLOBCM  is  required  to  couple  the  analysis  sub- 
routine directly  to  the  COPES/CONMIN.   All  variables,  which 
are  used  as  objective  function,  constraints  and  design  varia- 
bles, must  be  listed  in  the  common  statement  and  the  statement 
must  appear  in  each  subroutine  the  variables  are  used  in. 
It  is  used  by  the  optimizer  as  a  catalog  to  identify  where 
the  design  variables,  objective  function  and  constraints  are, 
and  what  purpose  they  fulfill. 
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In  order  to  execute  the  COPES/CONMIN  program  it  is 
necessary  to  provide  formatted  data  for  COPES/CONMIN,  followed 
by  data  for  the  ANALIZ  program.   Further  details  for  the 
COPES  data  are  explained  in  [Ref .  1] . 

D.   DESIGN  PAFAMETERS 

The  design  variables,  constraints  and  objective  function 
used  in  the  deisign  process  are  discussed  in  this  section. 

1 .  Design  Variables 

The  deisign  variables  are  plate  thicknesses,  t.,  and 
stiff ener  dimtmsions,  being  the  web  heights  h.,  the  web 
thicknesses  t,,.,  the  flange  widths  W.,  and  the  flange  thick- 
nesses tj..,    i  =  l,ns,  where  ns  is  the  number  of  the  different 
stiff  ener  dim(2nsion  sets.   Thus,  the  total  number  of  design 
variables  depends  on  the  sets  of  plate  thicknesses  and  the 
sets  of  different  stiff ener  dimensions. 

2 .  Objective  Function 

Total  structure  volume  is  considered  as  the  objec- 
tive function  in  the  design  process. 


m         ns 


Volume   =    y   V  .  +   y   N.  V  . 
i=l   P^    i=l   ^   ^^ 


where  V  .  is  the  volume  of  the  ith  plate  element,  m  is  the 
pi  ^ 

number  of  plate  elements,  V  .  is  the  volume  of  the  jth  stiff ener 
set,  and  N. is  the  number  of  stiff eners  in  this.   The  number 
of  stiff ener  sets  is  ns . 
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3.   Constraints 

Design  constraints  are  Von  Mises  maximum  stress, 
nodal  displacements,  height  to  thickness  ratios  of  frame 
webS/  and  width  to  thickness  ratios  of  frame  flanges. 

Stress: 


a../cr     -1<   0        i  =  l,ne 

il   max      —  .     '  . 

->  3  =  l^n: 

a  ■  ./o  -1<   0        i=l,ne 

11   max      —  .     '  . 

u       I 
where  a  ■  ■    and  a . .  are  respectively  the  upper  and  lower  sur- 
i:       i:        f  ji 

face  stress  at  the  node  j  of  the  element  i,  and  a     is  the 

-'  max 

maximum  allowable  stress.   The  number  of  elements  is  ne  and 
nj  is  the  number  of  nodal  points. 
Displacements: 


u.  ./u    -  1   <   0        i  =  l,nj 
^^   max      -  .  ^  ^^^^ 


where  u. .  is  the  displacement  at  node  i  in  the  coordinate 

direction  j  and  u    is  the  maximum  allowable  displacement 
-"      max  ^ 

in  coordinate  direction  j .   The  number  of  coordinates  is 
no. 

Height  to  thickness  ratio  of  frame  webs: 

VH^i  -  5   1  0  i  =  l,ns 

1  -  h./t,,.   <   0 
1   Wi   — 
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where  h.  and  t^.  are  respectively  the  web  height  and  web 
thickness  of  ith  stiffener  set. 

Width  to  thickness  ratios  of  flanges: 


w./t^.  -  5   <   0  i  =  l,ns 

1   fi      — 


1  -  w./t..   <   0 
i'^  fi   — 


where  w.  and  t^ .  are  respectively  the  width  and  the  thickness 
1      fi        f  J 

of  flange  of  the  ith  stiffener  set. 
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IV.   DESIGN  EXAMPLES 

This  chapter  will  present  design  examples  of  unstiffened 
and  stiffened  plates.   All  calculations  were  carried  out  on 
the  IBM  3033  model  370  computer  at  the  Naval  Postgraduate 
School.   Material  constants  used  for  all  examples  are  Young's 
modulus,  E  =  3  >:  10   psi,  and  Poisson's  ratio,  \>   =    0.3.   Due 
to  symmetry  only  a  quarter  of  the  plate  was  modelled  in  each 
example . 

A.   DESIGN  CASE  I — THE  UNSTIFFENED  PLATE 

A  quarter  of  the  unstiffened  plate  simply  supported,  as 
shown  in  Fig.  4.1,  was  modelled  using  a  4  x  4  mesh  of  the 
plate  elements.   There  are  10  design  variables  as  illustrated 
in  Fig.  4.2.   The  volume  of  the  plate  was  minimized  subject 
to  a  total  of  378  constraints  being  stress  and  nodal  deflec- 
tion constraints,  with  a    =  20,000  psi  and  u    =  1/2  t. 

max  ^        max 

1 .   Case  lA  with  A  Concentrated  Load 

A  concentrated  load  of  1,000  lbs  was  applied  at  the 
center  of  the  whole  plate.   The  final  optimization  results 
show  that  there  are  3  critical  deflection  constraints  of 
the  nodes  3,  4  and  10,  and  that  the  plate  thicknesses  of  the 

diagonal  elements  are  much  thicker  than  those  of  off-diagonal 

3 
elements.   The  final  volume  is  363.916  in  .   The  results  of 

this  case  are  summarized  in  Table  I  for  comparison  with  the 

other  cases. 
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t  :  plate  thickness 
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Figure  4.1    Design  Case  I — The  Unstiffened  Plate  Model 
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0)  ! 


Figure  4.2    Design  Case  I — Design  Variables 
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TABLE  I 
Summary  of  Design  Case  I 


Variable 

t(l) 

t(2) 

t(3) 

t(4) 

t(5) 

t(6) 

t(7) 

t(8) 

t(9) 

t(10) 

Vol 


Concentrated 

Load 

1. 

,  302 

1. 

.0  31 

0, 

.573 

1. 

.059 

0. 

.120 

0, 

.120 

0. 

.120 

0, 

.120 

0. 

.120 

0, 

.120 

363, 

.916 

Distributed 
Load 

0. 

.352 

0. 

.263 

0. 

.589 

0. 

.799 

0. 

.086 

0. 

.085 

0. 

.086 

0, 

.086 

0. 

.086 

0. 

.085 

170. 

.396 

t:   plata  thickness  (in) 
Vol:   volume  (in^) 


2 .   Case  IB  with  Uniform  Distributed  Load 

2 

A  uniform  distributed  load,  0.278  lbs/in   was  applied 

over  the  whole  plate.   The  final  optimization  results  show 
that  there,  are  6  critical  deflection  constraints  of  the 
nodes  3,  4,  5,  9,  10  and  15,  and  that  the  plate  thickness 
has  a  similar  trend  to  case  lA  where  diagonal  elements  are 

thicker  than  those  of  the  off -diagonal  elements.   The  final 

3 
volume  is  170.396  in  .   The  results  are  summarized  in  Table  I 

with  those  of  case  lA. 

B.   DESIGN  CASE  II — THE  STIFFENED  PLATE 

A  quarter  of  the  stiffened  plate  simply  supported,  as 
shown  in  Fig.  4.3,  was  modelled  4x4  mesh  of  the  plate  ele- 
ments and  12  stiff ener  elements. 
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Figure  4.3   Design  Case  II — The  Stiffened  Plate  Model 
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1 .   Case  IIA — Rectangular  Type  Frame  Stiffener 

There  are  17  design  variables:  a  plate  thickness, 
and  8  sets  of  stiffener  heights  and  thicknesses,  cis  tabu- 
lated in  Table  II.  The  volume  of  the  stiffened  p].ate  was 
minimized  subject  to  a  total  of  222  constraints:  stress, 
stiffener  height  to  thickness  ratio,  and  nodal  deflection 
constraints. 

a.  Case  IIAl  with  A  Concentrated  Load 

A  concentrated  load  of  1,000  lbs  v/as  cipplied  at 
the  center  of  the  whole  plate.   There  are  14  crit:.cal  con- 
straints:  the  height  to  thickness  ratios,  and  6  deflections 
of  the  nodes  4,  5,  10,  29,  30  and  31.   The  final  volume  is 
102.017  in  with  a  uniform  plate  thickness  of  0.0r)l  in.   The 
results  are  summarized  in  Table  III  for  comparing  with  the 
other  cases. 

b.  Case  IIA2  with  Uniform  Distributed  Load 

A  uniform  distributed  load,  0,278  lb/ in   was 
applied  over  the  whole  plate.   There  are  12  critical  con- 
straints:  the  height  to  thickness  ratios,  and  4  deflections 

of  the  plate  nodes  7,  8,  14  and  19.   The  final  volume  is 

3 
65.9  84  in   with  a  uniform  plate  thickness  of  0.059  in.   The 

results  of  this  case  are  summarized  in  Table  III  with  those 

of  the  case  IIAl.   The  volumes  of  the  above  two  cases  are 

significantly  reduced  in  comparison  to  those  of  the  unstif- 

fened  plate. 
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TABLE  II 
Design  Case  IIA — Design  Variables  and  Nodal  Connectivity 


element 
number 

design 
variables 

Nl 

node 
N2 

number 
N3 

N4 

Pl 

t 

1 

6 

7 

2 

p2 

t 

2 

7 

8 

3 

P3 

t 

3 

8 

9 

4 

p4 

t 

4 

9 

10 

5 

P5 

t 

6 

11 

12 

7 

p6 

t 

7 

12 

13 

8 

P7 

t 

8 

13 

14 

9 

p8 

t 

9 

14 

15 

10 

p9 

t 

11 

16 

17 

12 

plO 

t 

12 

17 

18 

13 

pll 

t 

13 

18 

19 

14 

pl2 

t 

14 

19 

20 

15 

pl3 

t 

16 

21 

22 

17 

pl4 

t 

17 

22 

23 

18 

pl5 

t 

18 

23 

24 

19 

pl6 

t 

19 

24 

25 

20 

si 

^^ 

r 

,h(l) 

30 

35 

s2 

\: 

.2] 

,h(2) 

35 

36 

s3 

Hi^ 

3' 

,h(3) 

36 

37 

s4 

Ht< 

A 

,h(4) 

37 

38 

s5 

\' 

[5 

,h(5) 

29 

30 

s6 

S;< 

.6] 

,h(6) 

28 

29 

s7 

"^ 

J' 

,h(7) 

27 

28 

s8 

^ 

.8; 

,h(8) 

26 

27 

s9 

^ 

;5; 

>,h(5) 

30 

31 

slO 

S^ 

;6; 

,h(6) 

31 

32 

sll 

S^ 

:7^ 

,h(7) 

32 

33 

sl2 

^ 

:8 

,h(8) 

33 

34 
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TABLE  III 
Summary  of  Design  Case  IIA 


variable 


^ 

a) 

s^ 

[2) 

^w' 

.3) 

S7< 

.4) 

^W< 

5) 

W 

.6) 

S.< 

7) 

S^< 

8) 

h  ( 

1) 

h  ( 

2) 

h  ( 

3) 

h  ( 

4) 

h  ( 

5) 

h  ( 

6) 

h  ( 

7) 

h 

.8) 

Vol 


concentrated 

load 

0 

.051 

0. 

.438 

0 

.318 

0. 

553 

0 

673 

0 

097 

0. 

096 

0. 

095 

0, 

095 

2, 

189 

1. 

563 

2. 

753 

3 

343 

0. 

486 

0 

.482 

0 

.478 

0. 

477 

102 

.017 

distributed 
load 

0 

.059 

0 

.209 

0 

.168 

0 

.267 

0 

.344 

0 

.138 

0 

.164 

0 

171 

0 

.155 

1 

.044 

0 

.832 

1 

.329 

1 

.709 

0 

.687 

0 

.821 

0 

.851 

0 

.776 

65 

.984 

t 

h 
Vol 


plate  thickness  (in) 
frame  thickness  (in) 
frame  height  (in) 
volume  (in  ) 
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2.   Case  IIB — T  Type  Frame  Stiffener 

There  are  33  design  variables:  a  plate  thickness,  8 
sets  of  the  web  heights  and  thicknesses,  and  8  sets  of  the 
flange  widths  and  thicknesses,  as  tabulated  in  Table  IV. 
The  volume  of  the  stiffened  plate  was  minimized  subject  to 
a  total  of  230  constraints:   stress,  web  height  to  thickness 
ratio,  flange  width  to  thickness  ratio,  and  nodal  deflection 
constraints. 

a.  Case  IIBl  with  A  Concentrated  Load 

A  concentrated  load  of  1,000  lbs  was  applied  at 
the  center  of  the  whole  plate.   There  are  12  critical  con- 
straints:  the  web  height  to  thickness  ratios,  the  3rd  and 
8th  flange  width  to  thickness  ratios,  and  deflection  con- 
straints at  the  node  5  of  the  element  p5  and  node  30  of  the 
element  si.   The  final  volume  is  55.091  in  with  a  uniform 
plate  thickness  of  0.023  in.   The  results  of  this  case  are 
summarized  in  Table  V  for  comparison  with  the  other  cases. 

b.  Case  IIB2  with  Uniform  Distributed  Load 

2 

A  uniform  distributed  load,  0.278  lbs/in  was 

applied  over  the  whole  plate.   There  are  16  critical  con- 
straints:  the  web  height  to  thickness  ratios  except  the 
8th  one,  the  3rd  through  8th  flange  width  to  thickness  ratios 

except  the  6th  one,  and  4  deflections  of  the  plate  nodes  7, 

3 
8,  14  and  19.   The  final  volume  is  64.414  in   with  a  uniform 

plate  thickness  of  0.06  0  in.   This  volume  is  greater  than  the 

volume  of  the  above  concentrated  load  case,  being  different 

from  the  other  cases  which  the  latter  is  greater  than  the 
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TABLE  IV 
Design  Case  IIB — Design  Variables  and  Nodal  Connectivity 


element 
number 

Pl 
P2 
p3 
p4 

P5 
p6 

P7 

p8 

p9 

plO 

pll 

pl2 

pl3 

pl4 

pl5 

pl6 

si 

s2 

s3 

s4 

s5 

s6 

s7 

s8 

s9 

slO 

sll 

sl2 


^w^i-^ 


2 
■3 

:4 
5: 

6 

7: 

8 

5: 
6; 

7; 

8 


design 
variables 


t 

t 

t 

t 

t 

t 

t 

t 

t 

t 

t 

t 

t 

t 

t 

t 

h(l 

h(2 

h(3 

h(4 

h(5 

h(6 

h(7 

h(8 

h(5 

h(6 

h(7 

h(8 


,w(l 
,w(l 
,w{3 
,w(4 
,w(5 
,w(6 
,w(7 
,w(8 
,w(5 
,w(6 
,w(7 
,w(8 


node 

number 

Nl 

N2 

N3 

N4 

1 

6 

7 

2 

2 

7 

8 

3 

3 

8 

9 

4 

4 

9 

10 

5 

6 

11 

12 

7 

7 

12 

13 

8 

8 

13 

14 

9 

9 

14 

15 

10 

11 

16 

17 

12 

12 

17 

18 

13 

13 

18 

19 

14 

14 

19 

20 

15 

16 

21 

22 

17 

17 

22 

23 

18 

18 

23 

24 

19 

19 

24 

25 

20 

30 

35 

35 

36 

36 

37 

37 

38 

29 

30 

28 

29 

27 

28 

26 

27 

30 

31 

31 

32 

32 

33 

33 

34 
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TABLE  V 
Summary  of  Design  Case  IIB 


variable 

t 

tw< 

'D 

^W 

[2) 

^w' 

3) 

^w< 

4) 

^w< 

5) 

tw< 

.6) 

tw< 

7) 

^w< 

8) 

h  ( 

1) 

h  ( 

2) 

h  ( 

3) 

h  ( 

4) 

h  ( 

5) 

h 

.6) 

h 

:7) 

h 

.8) 

^f 

(1) 

tf 

(2) 

tf 

(3) 

4 

i4) 

tf 

[5) 

tf 

[6) 

^f 

(7) 

t^ 

(8) 

concentrated 
load 

0. 

023 

0. 

133 

0, 

077 

0. 

137 

0. 

162 

0. 

381 

0. 

355 

0. 

316 

0. 

274 

0. 

659 

0. 

386 

0. 

680 

0. 

.805 

1 

.902 

1 

767 

1 

.570 

1 

.  368 

0 

.021 

0 

.028 

0 

.011 

0 

.026 

0 

.876 

0 

.799 

0 

.671 

0 

.460 

distributed 
load 

0. 

060 

0. 

149 

0. 

126 

0. 

197 

0. 

242 

0. 

103 

0. 

114 

0. 

115 

0. 

122 

0. 

736 

0 

.631 

0 

.976 

1 

.201 

0 

.509 

0 

.566 

0 

.569 

0 

.602 

0 

.262 

0 

.183 

0 

.308 

0 

.460 

0 

.209 

0 

.231 

0 

.257 

0 

.201 
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TABLE  V  (Cont'd) 


variable 

w 

:i) 

w  1 

2) 

w  < 

.3) 

w  ( 

4) 

w  ( 

5) 

w  ( 

6) 

w  ( 

7) 

w  ( 

8) 

Vol 


concentrated 

load 

0 

.083 

0. 

054 

0 

.057 

0. 

074 

0 

.891 

0. 

814 

0 

685 

0 

.464 

55 

091 

t:  place  thickness  (in) 

t.,:  web  thickness  (in) 
w 

h:  web  height  (in) 

t^:  flange  thickness  (in) 

w:  flange  width  (in) 

Vol:  volume  (in-^) 


distributed 
load 

0 

.2  78 

0 

.18  7 

0 

.307 

0 

.4  60 

0 

.210 

0 

.2  35 

0 

.258 

0 

.201 

64 

.414 
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former.   This  suggests  that  the  stiff ener  type  and  arrange- 
ment used  here  is  particularly  efficient  for  the  concentrated 
load  case.   The  results  of  this  case  are  summarized  in  Table 
V  with  those  of  the  case  IIBl.   The  volumes  of  the  above 
cases  are  significantly  reduced  with  comparison  to  those  of 
unstiffened  plates.   The  volume  of  T  type  stiff ener  case 
is  slightly  smaller  than  that  of  rectangular  case  in  the 
distributed  load.   On  the  other  hand,  the  volume  of  T  type 
stiffener  case  is  much  smaller  than  the  rectangular  case 
for  the  concentrated  load. 

C.   DESIGN  CASE  III — THE  HATCH  COVER 

A  quarter  of  a  hatch  cover  was  modelled  4x4  mesh  of  the 
plate  elements  and  24  stiffener  elements,  as  shown  in  Fig. 
4.4. 

1.   Case  IIIA-- Rectangular  Type  Frame  Stiffener 

There  are  9  design  variables:   a  plate  thickness, 
and  4  sets  of  stiffener  heights  and  thicknesses,  as  tabulated 
in  Table  VI.   The  volume  of  the  hatch  cover  was  minimized 
subject  to  a  total  of  2  74  constraints:   stress,  height  to 
thickness  ratio,  and  nodal  deflection  constraints. 

a.   Case  IIIAl  with  A  Concentrated  Load 

A  concentrated  load  of  1,000  lbs  was  applied  at 
the  center  of  the  whole  plate.   There  are  4  critical  con- 
straints:  the  2nd  and  4th  height  to  thickness  ratios,  and 
the  upper  and  lower  stresses  of  the  node  5  of  the  element 
p4 .   The  dimensions  of  the  2nd  stiffener  set  are  negligibly 
small  with  comparison  to  those  of  the  other  stiffener  sets. 
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Figure  4.4    Design  Case  III — The  Hatch  Cover  Model 
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TABLE  VI 
Design  Case  IIIA — Design  Variables  and  Nodal  Connectivity 


element 

design 

number 

variables 

Pl 

t 

p2 

t 

P3 

t 

p4 

t 

P5 

t 

p6 

t 

P7 

t 

p8 

t 

p9 

t 

plO 

t 

pll 

t 

pl2 

t 

pl3 

t 

pl4 

t 

pl5 

t 

pl6 

t 

si 

l/2t^(3) ,h(3) 

s2 

l/2t^(3) ,h(3) 

s3 

t^(l) ,h(l) 

s4 

t^(l) ,h(l) 

s5 

tw(4),h(4) 

s6 

t^(2),h(2) 

s7 

t^(l),h(l) 

s8 

t^(4),h(4) 

s9 

t^(2),h(2) 

slO 

t^(l) ,h(l) 

sll 

Ht(3)  ,h(3) 

sl2 

t,,(3),h(3) 

node  number 


Nl 

N2 

N3 

N4 

1 

6 

7 

2 

2 

7 

8 

3 

3 

8 

9 

4 

4 

9 

10 

5 

6 

11 

12 

7 

7 

12 

13 

8 

8 

13 

14 

9 

9 

14 

15 

10 

11 

16 

17 

12 

12 

17 

18 

13 

13 

18 

19 

14 

14 

19 

20 

15 

16 

21 

22 

17 

17 

22 

23 

18 

18 

23 

24 

19 

19 

24 

25 

20 

26 

27 

27 

28 

28 

29 

29 

30 

26 

31 

28 

32 

30 

33 

31 

34 

32 

36 

33 

38 

34 

35 

35 

36 
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TABLE   VI     (Cont'd) 


element 

design 

node 

number 

number 

variableis 

Nl 

N2 

N3 

N4 

si  3 

t^,(2) 

,h{2) 

36 

37 

sl4 

S,(2). 

,lrA2) 

37 

38 

si  5 

tw(4) 

rh.(4) 

34 

39 

sl6 

t.,(3), 

h(3) 

36 

40 

sl7 

l/2t^,(3)< 

.h(3) 

38 

41 

si  8 

^-w(^) 

,h(4) 

39 

42 

si  9 

\.(3), 

-h(3) 

40 

44 

s20 

l/2t^,(3) 

,h(3) 

41 

46 

s21 

t,,(4) 

rh(4) 

42 

43 

s22 

tw(4) 

.h(4) 

43 

44 

s23 

H^(^)' 

-h(4) 

44 

45 

s24 

t,,(4). 

-h(4) 

45 

46 
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This  suggests  that  the  1st  stiffener  set  is  stiff  enough  to 
compensate  the  elements  attached  to  the  2nd  stiffener  set. 
The  final  volume  is  323.639  in  with  a  uniform  plate  thick- 
ness of  0.359  in.   The  results  of  this  case  are  summarized 
in  Table  VII. 


variable 


t^(l) 

t^(4) 
h  (1) 
h  (2) 
h  (3) 
h  (4) 

Vol 


W 
h 

Vol 


TABLE  VII 
Summary  of  Design  Case  IIIA 

concentrated 
load 

0,359 

0.811 

0.001 

0.084 

0.015 

1.306 

0.006 

0.162 

0.073 

323.639 

plate  thickness  (in) 
frame  thickness  (in) 
frame  height  (in) 
volume  (in  ) 


distributed 
load 

0.059 

0.863 

0.002 

0.518 

0.009 

1.750 

0.011 

1.484 

0.019 

87.314 


b.   Case  IIIA2  with  Uniform  Distributed  Load 

2 

A  iiniform  distributed  load,    0.278    lbs/in     was 

applied  over  the  whole  plate.   There  are  6  critical  constraints 
the  2nd  height  to  thickness  ratio,  and  5  deflections  of  the 
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plate  nodes  8,  9  and  14  and  of  stiff ener  nodes  32  and  37. 
With  the  same  suggestion  as  the  concentrated  load  case  the 

dimensions  of  the  2nd  stiff ener  set  are  negligibly  small. 

3 
The  final  volume  is  8  7.314  in  with  a  uniform  plate  thickness 

of  0.059  in.   The  design  results  of  the  above  both  cases 

are  summarized  in  Table  VII  with  those  of  case  IIIAi.   The 

volumes  of  the  above  two  cases  are  significantly  reduced 

with  comparison  to  those  of  unstiffened  plates. 

2.   Case  IIIB — T  Type  Frame  Stiffener 

There  are  17  design  variables:   a  plate  thickness,  4 
sets  of  the  web  heights  and  thicknesses,  and  4  sets  of  the 
flange  widths  and  thicknesses,  as  tabulated  in  Table  VIII. 
The  volume  of  the  hatch  cover  was  minimized  subject  to  a  t^tal 
of  2  78  constraints:   stress,  web  height  to  thickness  r?.tio, 
flange  width  to  thickness  ratio,  and  nodal  deflection 
constraints. 

a.   Case  IIIBl  with  A  Concentrated  Load 

A  concentrated  load  of  1,000  lbs  was  appl.^ed  at 
the  center  of  the  whole  plate.   There  are  3  critical  con- 
straints:  the  3rd  flange  width  to  thickness  ratios,  and  2 
deflections  of  the  plate  node  5  and  stiffener  node  30. 
Negligibly  small  are  the  whole  dimensions  of  the  2nd  stiffener 
set  and  the  flange  dimensions  of  the  3rd  stiffener  set  with 

comparison  to  the  other  dimensions.   The  final  volume  is 

3 
95.883  in  with  a  uniform  plate  thickness,  0.083  in. 
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TABLE  VIII 
Design  Case  IIIB — Design  Variables  and  Nodal  Connectivity 


element 

di 

2  sign 

node 

number 

nioinber 

van 

.ables 

Nl 

N2 

N3 

N4 

Pl 

t 

1 

6 

7 

2 

p2 

t 

2 

7 

8 

3 

P3 

t 

3 

8 

9 

4 

p4 

t 

4 

9 

10 

5 

P5 

t 

6 

11 

12 

-7 

p6 

t 

7 

12 

13 

8 

P7 

t 

8 

13 

14 

9 

p8 

t 

9 

14 

15 

10 

p9 

t 

11 

16 

17 

12 

plO 

t 

12 

17 

18 

13 

pll 

t 

13 

18 

19 

14 

pl2 

t 

14 

19 

20 

15 

pl3 

t 

16 

21 

22 

17 

pl4 

t 

17 

22 

23 

18 

pl5 

t 

18 

23 

24 

19 

pl6 

t 

19 

24 

25 

20 

si 

l/2t^< 

3) 

,h(3^  , 

^f 

[3] 

,l/2w(3) 

26 

27 

s2 

l/2t^ 

;3) 

,h(3)  , 

[3] 

,l/2w(3) 

27 

28 

s3 

^W 

[1) 

,h{l}  , 

"-f 

(i; 

,w(l) 

28 

29 

s4 

^W 

[1) 

,h(l)  , 

■1- 

'  -f 

(i; 

,w(l) 

29 

30 

s5 

^W 

[4) 

,h{4) , 

'tf 

(4 

,w(4) 

26 

31 

s6 

^ 

(2) 

rh{2) , 

[2; 

,w(2) 

28 

32 

s7 

^W 

(1) 

,h(l) . 

'^f 

(1 

,w(l) 

30 

33 

s8 

^W 

(4) 

,h(4)  , 

'^f 

(4 

,w(4) 

31 

34 

s9 

^W 

[2) 

,h(2)  , 

'tf 

(2 

,w(2) 

32 

36 

slO 

^W 

(1) 

,h(l) 

'tf 

[1 

l,w(l) 

33 

38 

sll 

^W 

(3) 

,h(3) 

'tf 

(3^ 

,w(3) 

34 

35 

sl2 

tT.T 

(3) 

,h(3) 

'^f 

(3^ 

,w(3) 

35 

36 
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TABLE  VIII  (Cont'd) 


element 
number 

sl3 
sl4 
sl5 
sl6 
sl7 
sl8 
sl9 
s20 
s21 
s22 
s23 
s24 


design 
variables 


l/2t 


1/2 1 


.„(2), 

.h(2) 

,tj( 

2), 

W(2) 

.h(2), 

tf 

:2), 

.„(4) 

,h(4) , 

■*^f 

,4)  , 

.„(3) 

,h(3)  , 

tf 

3)  , 

„(3), 

h(3)  , 

tf 

;3) , 

•„(4) 

.h(4), 

tf 

;4), 

„(3), 

h(3)  , 

tf 

;3), 

•„(3), 

.h(3)  , 

■tf 

.3)  , 

^(4), 

.h(4), 

■tf' 

,4)  , 

„(4), 

h(4), 

tf 

4), 

.„(4), 

h(4)  , 

■tf' 

4)  , 

.„(4), 

,h(4), 

tf 

14), 

w(2) 
w(2) 
w(4) 
w(3) 


Nl 

36 
37 
34 
36 


l/2w(3)  38 

w(4)  39 

w(3)  40 

l/2w(3)  41 

w(4)  42 

w(4)  43 

w(4)  44 

w(4)  45 


node  number 
N2    N3    N4 

37 
38 
39 

40 
41 
42 
44 
46 
43 
44 
45 
46 
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b.   Case  IIIB2  with  Uniform  Distributed  Load 

2 
A  uniform  distributed  load,  0.2  78  Ibs/m   was 

applied  over  the  whole  plate.   There  are  9  critical  con- 
straints:  the  3rd  web  height  to  thickness  ratio,  the  3rd 
flange  width  to  thickness  ratio,  4  deflections  of  the  plate 
nodes  8,  9,  13  and  14,  and  3  deflections  of  the  stiff ener 
nodes  32,  36  and  37.   Negligibly  small  are  the  v/hole  dimen- 
sions of  the  2nd  stiff ener  set,  and  the  flange  dimensions  o: 
the  3rd  and  4th  stiffener  sets  with  comparison  to  the  ether 
dimensions.   This  suggests  that  the  1st  stiffener  set  is 
stiff  enough  to  compensate  the  elements  attached  to  the  2nd 
stiffener  set.   Additionally  it  can  be  suggested  that  the 
3rd  and  4th  stiffener  sets  are  stiff  enough  to  support  the 

given  load  with  only  the  rectangular  type  ones.   The  final 

3 
volume  is  76.4  36  in   with  a  uniform  plate  thickness,  0,0  62 

in.   The  results  of  the  above  both  cases  are  sumonarized  in 
Table  IX.   The  volumes  of  the  above  two  cases  are  signifi- 
cantly reduced  with  comparison  to  those  of  the  unstiffened 
plates  and  the  volume  of  T  type  stiffener  case  is  much 
smaller  than  that  of  the  rectangular  type  case. 
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TABLE  IX 
Summary  of  Design  Case  IIIB 

variable  concentrated  distributed 

load  load 

t  0.083  0.062 

0.284  0.284 

0.003  0.003 

0.323  0.318 

0.020  0.020 

1.201  1.200 

0.010  0.010 

1.433  1.430 

0.060  0.060 

1.001  1.000 

0.003  0.003 

0.000  0.001 

0.010  0.010 

1.110  1.110 

0.010  0.010 

0.002  0.003 

0.030  0.001 

Vol  95.883  76.436 

plate  thickness  (in) 
web  thickness  (in) 
web  height  (in) 
flange  thickness  (in) 
flange  width  (in) 
volume  (in  ) 


S*' 

1) 

tw 

[2) 

'w 

.3) 

*w 

:4) 

h    ( 

.1) 

h    ( 

2) 

h    ( 

3) 

h    ( 

4) 

tf' 

1) 

•^f' 

2) 

tf 

'3) 

tf' 

4) 

w 

:i) 

w 

.2) 

w 

[3) 

w 

(4) 

t: 

S«  = 

h: 

tf  = 

w: 

Vol: 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

A.   CONCLUSIONS 

A  finite  element  analysis  program  was  developed  and 
coupled  to  an  optimizer  for  stiffened  plate  design.   Stif- 
fened plates  were  designed  for  minimum  volume  with  a  concen- 
trated load. 

Five  important  conclusions  are  made  from  considering 
several  design  examples. 

First,  a  stiffened  plate  is  much  more  efficient  than 
the  unstiffened  plate  for  the  optimum  volume  design. 

Second,  it  is  required  to  attach  the  larger  stiff eners 
to  the  diagonal  elements  of  the  plate  as  compared  to  the 
off -diagonal  elements. 

Third,  the  T  type  stiffeners  provide  a  lower  volume 
design  than  the  rectangular  type  stiffeners. 

Fourth,  the  stiffened  plate  with  the  diagonally  attached 
stiffeners  are  much  more  efficient  than  those  with  the 
rectangularly  attached  stiffeners. 

Fifth,  the  stiffened  plate  attached  diagonally  with  the 
T  type  stiffeners  is  particularly  efficient  for  the  concen- 
trated load  case. 

These  conclusions  are  based  on  the  design  of  simply 
supported  plates  and  the  conclusions  could  be  quite  differ- 
ent if  other  boundary  conditions  are  considered. 


5  7 


B.   RECOMMENDATIONS 

The  following  recommendations  may  be  useful  for  future 
work. 

1.  The  analysis  program  should  be  modified  to  a 
generalized  program  with  a  variety  of  elements. 

2.  The  analysis  program  should  be  extended  to 
consider  the  dynamic  problem. 

3.  The  analysis  program  should  be  extended  to 
consider  the  multiple  loading  conditions. 

4 .  Routines  should  be  added  to  calculate  gradients 
analytically. 

5.  Buckling  constraints  of  the  plate  and  stiff ener 
elements  should  be  added  in  the  optimization 
process . 

6 .  Frequency  constraints  should  be  added  in  the 
optimization  process. 

7.  The  analysis  program  shoitild  be  modified  to  reduce 
the  computing  time  during  the  optimization  process 


58 


LIST  OF  REFERENCES 


1.  Madsen,  L.E.,  and  Vanderplaats ,  G.N.,  COPES — A  Fortran 
Control  Program  for  Engineering  Synthesis ^  NPS  69-81- 
003,  Naval  Postgraduate  School,  Monterey,  CA. ,  March 
1982. 

2.  Vanderplaats,  G.N,,  CONMIN — A  Fortran  Program  for 
Constrained  Function  Minimization--User * s  Manual, 
NASA  TMX-6  2,2  82,  August  19  73. 

3.  Przemieniecki,  J.S.,  Theory  of  Matrix  Structural  Analysis, 
McGraw-Hill,  1968,  pp.  96-125. 

4.  Miller,  R.E.,  "Reduction  of  the  Error  in  Eccentric  Beam 
Modelling, "  International  Journal  for  Numerical  Methods 
in  Engineering,  Vol.  15,  1980,  pp.  575-582. 

5.  Fox,  R.L.,  Optimization  Method  for  Engineering  Design, 
Addison-Wesley  Publishing  Co.,  Reading,  Mass.,  1971. 

6.  Himmelblau,  D.M.,  Applied  Nonlinear  Programming,  McGraw- 
Hill  Book  Co.,  New  York,  19  72. 

7.  Fletcher,  D.  and  Reeves,  C.K.,  "Function  Minimization 
by  Conjugate  Gradients,"  Computer  Journal,  Vol.  7, 
No.  2,  July  1964,  pp.  149-154. 

8.  Zoutendijk,  G.,  Methods  of  Feasible  Directions, 
Elsevier  Publishing  Co.,  Amsterdam,  1960. 

9.  Vanderplaats,  G.N.,  and  Moses,  F. ,  "Structural  Optimiza- 
tion by  Methods  of  Feasible  Directions,"  Computers  and 
Structures ,  Vol.  3,  July  1973. 

10.  Vanderplaats,  G.N.,  Hicks,  Raymond  H.,  and  Murman,  E.M. , 
"Application  of  Numerical  Optimization  Techniques  to 
Airfoil  Design,"  Aerodynamic  Analysis  Requiring  Advanced 
Computers,  NASA  SP-347,  Part  II,  March  1975,  pp.  749- 
768. 

11,  Vanderplaats,  G.N.,  Automated  Optimization  Techniques 
for  Aircraft  Synthesis,  AIAA  Paper  No.  76-909,  AIAA 
Aircraft  Systems  and  Technology  Meeting,  Dallas,  Texas, 
September  19  76. 


59 


INITIAL  DISTRIBUTION  LIST 


No.  Copies 


1.  Defense  Technical  Information  Center  2 
Cameron  Station 

Alexandria,  Virginia   22314 

2.  Library,  Code  014  2  2 
Naval  Postgraduate  School 

Monterey,  California   9  394  0 

3.  Library  2 
Republic  of  Korea  Naval  Academy 

Chinhae,  Republic  of  Korea 

4.  Department  Chairman,  Code  69  1 
Department  of  Mechanical  Engineering 

Naval  Postgraduate  School 
Monterey,  California   93940 

5.  Dr.  Garret  N.  Vanderplaats  10 
Code  69Vn 

Department  of  Mechanical  Engineering 
Monterey,  California   93940 

().   Dr.  Young  S.  Shin  1 

Code  69Sg 

Department  of  Mechanical  Engineering 
Monterey,  California   9  3940 

7.   LCDR  Chil  Sung  Park,  ROKN  5 

Personnel  Education  Division 
Naval  Headquarters 
Seoul,  Republic  of  KOREA 


60 


Thesis 
PI.5737 
c.l 


u 


'2'* 


Park 


Optimum  design 
stiffened  plater 
the  finite  ^^"^^ 

method 


lhesP15737 

Optimum  design  of  stiffened  plates 


usmg 


3  2768  001  97195  5 

DUDLEY  KNOX  LIBRARY 


